Testing model independent modified gravity with future large scale surveys 
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Model-independent parametrisations of modified gravity have attracted a lot of attention over 
the past few years and numerous combinations of experiments and observables have been suggested 
to constrain the parameters used in these models. Galaxy clusters have been mentioned, but not 
looked at as extensively in the literature as some other probes. Here we look at adding galaxy 
clusters into the mix of observables and examine how they could improve the constraints on the 
modified gravity parameters. In particular, we forecast the constraints from combining Planck 
satellite Cosmic Microwave Background (CMB) measurements and Sunyaev-Zeldovich (SZ) cluster 
catalogue with a DES-like Weak Lensing (WL) survey. We find that cluster counts significantly 
improve the constraints over those derived using CMB and WL. We then look at surveys further 
into the future, to see how much better it may be feasible to make the constraints. 

PACS numbers: 04.80.Cc,98.65.Cw,98.80.-k 



I. INTRODUCTION 

Einstein's General Relativity (GR) is one of the princi- 
pal ingredients of modern cosmology. Indeed, it could be 
argued that it was only with the development of GR that 
cosmology really became a part of physics. Nonetheless, 
it is our job as physicists to continue to test even the 
most fundamental pillars of cosmology in order to refine, 
improve and further justify our model of the universe. 
There are also fundamental reasons for considering dif- 
ferent theories of gravity: GR is inconsistent with quan- 
tum mechanics and the search for a theory of 'Quantum 
Gravity' is one of the holy grails of modern physics. 
It has proved to be difficult to test GR outside of the solar 
system, particularly as the effects of a different theory of 
gravity could be degenerate with behaviour induced by 
different constituents of the universe. This is the case 
with current observations that suggest the presence of 
some form of dark matter and dark energy. Dark matter 
is required to explain galaxy rotation curves, galaxy lens- 
ing, nucleosynthesis, acoustic oscillations in the CMB, 
and the growth of structure in the Universe to name a 
few. The requirement for dark energy is underpinned by 
observations of the background expansion rate and large 
scale structure measurement. 

The question has often been raised as to whether these 
effects could be due to a modified gravity theory. Propos- 
als for such a theory include a number of /(i?) theories, 
the Dvali Gabadadze Porrati (DGP) model [T], confor- 
mal gravity [2], Modified Newtonian dynamics (MOND) 
[3] and its covariant, relativistic extensions [1], and Ein- 
stein Aether theories [S], (see also the recent review [5j 
for an exhaustive list of candidates). It is possible, if GR 
is indeed the correct theory of gravity, that this debate 
will only be settled by the non-gravitational detection of 
the dark matter and/or dark energy. However, since we 
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do not know whether or not GR is indeed the correct 
theory of gravity, it seems reasonable to consider which 
observations could allow us to detect modified gravity. 

More importantly, are there ways to test deviations from 
GR in a model independent way? There are several ad- 
vantages to a model independent approach; some alter- 
natives to GR do exist but there is no complete theory 
of, for example, quantum gravity to draw on. Also, there 
are no 'stand-out' candidates that are universally consid- 
ered to be strong alternatives. Another advantage of a 
model independent approach is that the results do not 
rely on model selection techniques such as x^ per degree 
of freedom or other ways of choosing between compet- 
ing theories. Thus, they also do not rely on us having 
the correct theory to hand. A result that is inconsistent 
with the GR based, concordance cosmology will be un- 
ambiguous and therefore a strong motivator to develop 
alternative theories, as well as possibly giving us a clue 
as to the nature of these theories. 

There have been many recent studies looking at model 
independent tests of the dark energy /cold dark matter 
(ACDM) paradigm in GR. These can be roughly split 
into two categories: consistency checks of the ACDM as- 
sumption [7H5] and those that introduce new parameters 
to evaluate the level of deviation from GR (see [6] for a 
complete review of the literature). These parametrisa- 
tions and consistency checks have considered the major- 
ity of cosmological observations: weak lensing, the CMB, 
particularly through its Integrated Sachs Wolfe (ISW) ef- 
fect, Baryon Acoustic Oscillations (BAO), SNla Super- 
novae luminosity distance observations, cluster counts, 
and galaxy redshift and peculiar velocity surveys. The 
purpose is normally to constrain both the background 
expansion history and the perturbations, or growth of 
structure, around the background, with some observ- 
ables being sensitive to both. The constraints are then 
combined. This approach works well because the growth 
and expansion are determined by the same quantities un- 
der the concordance cosmology and this is the basis for 
the consistency checks. A parameter that can be used 



without modifying gravity is 7, a parametrisation of the 
growth index [lOl [11] : 
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Here, D(k,z) = S{k, z)/S{k, z = 00) where d is the mat- 
ter density contrast. For ACDM, 7 — 0.55 is a good fit to 
the growth. With the modified gravity parametrisations, 
the parameters often relate to the two gravitational po- 
tentials ^ and $ appearing in the perturbation to the 
Friedmann Roberston Walker (FRW) metric. Since dif- 
ferent observations depend on different combinations of 
these potentials, combining several experiments gives the 
best constraints. 

Constraints from current data have been examined (see 
[5J for a review) and the general conclusion appears to 
be that the concordance cosmology is consistent with all 
of the current data. However, the data available today 
does not have enough constraining power to rule out even 
relatively significant modifications to gravity and they 
are certainly not precise enough to distinguish between 
modified gravity theories. Work has also gone into fore- 
casting future constraints on a number of theories and/or 
parametrisations of the modifications to GR [T^HT^ . 
The consensus is that future surveys will greatly improve 
prospects with the expectation that a number of theories 
competing with GR will be ruled out. 
Now is a good time to consider these issues as we are in an 
era where the observational front is rapidly advancing in 
the field. Most, if not all, of the current and future mea- 
surements can be brought to bear on the issue of modified 
gravity. We have had pioneering ground and space-based 
CMB experiments over the last 20 years and now await 
the data from the Planck satellite [20] . Weak lensing sur- 
veys of cosmic shear have also matured into a precise ob- 
servational tool [21] and future, planned surveys promise 
to bring these measurements to the fore front of the data 
landscape with large scale surveys underway or in the 
development stage (e.g. PANSTARRS ^, DES [13], 
and LSST 24J). In addition large scale surveys of galaxy 
redshifts have already been carried out (SDSS [2S], and 
2DFGRS 26 ) with even larger and deeper ones targeting 
BAO measurements to come. 

So far, we have only mentioned modified gravity as an 
alternative to dark energy. The other alternative under 
consideration is whether the assumptions of homogene- 
ity and isotropy are justified and hence whether using 
the FRW metric itself is justified. It is also important 
to note that measurements of the background expansion 
alone cannot distinguish between dark energy or modi- 
fied gravity [37] and it has also been argued that mea- 
surements of perturbations may also suffer from this de- 
generacy if sufficiently complex models of perturbed dark 
energy are allowed [25]. 

In this work, we will only consider modified gravity and 
will look at how well some future experiments, particu- 
larly DES and CMB measurements and SZ cluster counts 
from Planck and its successors, will combine to constrain 



certain modified gravity parameters that characterise po- 
tential deviations from GR. 

In this paper we investigate how the combination of fu- 
ture observations of CMB, weak lensing (WL), and clus- 
ter counts (CC) are able to constrain the model indepen- 
dent parametrisation of modified gravity theories. We re- 
strict ourselves to the simplest form of modified gravity 
with at most a linear redshift dependence in the mod- 
ification and no scale dependence. These assumptions 
are fairly restrictive in terms of the physical mechanism 
that could underpin the modified phenomenology but 
provide a simple starting point for investigation into fu- 
ture constraints from different observables. In Section HIl 
we briefly review the model independent parametrisation 
of modified gravity used in this work. In Section |III| we 
describe the three observables used in our forecasts. In 
Section lrvl we review the Fisher matrix formalism used in 
our calculations and the experimental 'Stages' considered 
in our forecasts. Our results are presented in Section \V\ 
and we conclude with a discussion of the results in Sec- 
tion IVH 



II. PARAMETRISED MODIFIED GRAVITY 

There has been considerable discussion regarding the 
best ways to address and parametrise deviations from 
GR. Of course, it is always possible to calculate observ- 
ables in the universe given a particular modified grav- 
ity model, and this has been done for several models: 
/(i?), DGP, Einstein Aether, TeVeS and conformal grav- 
ity have all been tested against cosmological data. How- 
ever, it is useful to consider general kinds of deviations 
from GR that can then be used as a 'null-test' of the 
current paradigm. In other words, if there is significant 
evidence that the parameters we consider differ from their 
values under GR, then there is clearly a case that the cur- 
rent model is wrong that is not predicated on a particular 
modified gravity model. 

Numerous sets of 'modified gravity parameters' (MGPs) 
have been suggested in the literature, see e.g. [2S] for a 
partial translation table and [T^ for a discussion of the 
differences with some of the parametrisations. Most of 
the parametrisations are phenomenological modifications 
to the Einstein equations and typically involve a param- 
eter relating to the strength of gravity and a parameter 
relating the two scalar potentials in the metric. 
In this work we will use two parameters, 77 and /i, follow- 
ing [3D]. The first, 77, is the ratio of the two potentials, 
77 = ^/^P. This will be roughly equal to unity in GR 
unless any of the particle species has large anisotropic 
stress, this is not expected to be the case unless a signif- 
icant amount of dark matter is made up of massive neu- 
trinos in conflict with Large Scale Structure (LSS) data. 
The second, fi, is a modification of the poisson equation, 
and is essentially a time and space dependent Newton's 
constant. Fourier expanding the spatial dependence with 
wavenumbers k and assuming isotropy, the modification 



of the Poisson equation is as follows 

k^-^{a, k) = -47rGaV(a, fc)p(a)A(a, k) , (2) 

where, a is the FRW scale factor, G is Newton's constant, 
p is the background density of cold dark matter and A 
is the gauge invariant density contrast given by 
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where the cold dark matter density contrast is defined as 
S = 6p{a, k)/p{a), v is velocity of the dark matter and H 
is the Hubble parameter. 

The GR limit is recovered when the function p is constant 
and equal to unity. Note that we are using the potential 
^ in the Poisson equation since this determines the ac- 
celeration of non-relativistic particles. These MGPs have 
been included as parameters in the modified Einstein- 
Boltzmann code MGCAMB 30] based on the well-known 
CAMB 31 package for integrating photon and matter per- 
turbations. This parametrisation can be related directly 
to particular models through the definition of /3 and A 
parameters [271150] . 

Our potentials are defined as scalar perturbations of a 
flat, FRW metric 



500 = -[l + 2*(f,i)] , 
5„- = a^{t)[1^2^{x,t)]S,^ 



(4) 



where we have made the conformal Newtonian gauge 
choice to fix the remaining two scalar degrees of freedom 
in the perturbed metric. In (E), ^ is the Newtonian po- 
tential and is responsible for the acceleration of massive 
particles whereas $ is the curvature potential, which also 
contributes to the acceleration of relativistic particles. 
We will consider a number of scenarios for the phe- 
nomenological modification to the standard gravitational 
force. We will assume that GR is valid up to a specified 
redshift. This is motivated by the stringent conditions 
set on any modifications to the standard paradigm by 
Big Bang Nucleosynthesis (BBN) observations of light 
element abundances in the early universe and of CMB 
anisotropics at recombination. Beyond the redshift where 
modifications take over we will assume either that the p 
and rj parameters remain constant with a given value 
different from unity or that they follow a simple time- 
dependence. 

This simple treatment raises two issues. Firstly, the time- 
dependence of the modification introduced in ([2| is not 
motivated by any dynamical considerations and may not 
be self-consistent within the framework of the underlying 
theory responsible for the departure from GR. This is an 
unavoidable problem in the phenomenological approach 
taken here ^ . Secondly our prescription does not allow 



^ However see |32| for recent work in defining self-consistent MGPs 
for expansions around GR. 



for any spatial dependence of the modifications. This 
can be justified in part by requiring the simplest modi- 
fications however all mechanisms which generate devia- 
tions from GR will have to include some form of screening 
or cut-off mechanism, at the very least, acting on scales 
close to solar system and below in order to satisfy labora- 
tory and solar system scale tests of GR. Indeed, all mech- 
anisms proposed so far have have some form of explicit 
screening mechanism acting either universal or environ- 
ment dependent length scales. For the purpose of this 
exercise we shall omit any screening mechanism, effec- 
tively assuming it acts on scales smaller than any probed 
by the observations included in this work. However if the 
screening scale is effective at even moderate comoving 
scales of a few Mpc then we would expect some sensitiv- 
ity of our results to this assumption, particularly where 
we have included cluster counts. We leave for future work 
the determination of the effect of scale dependent mod- 
ifications on cluster count predictions which necessarily 
require the study of modified N-body simulations [55] . 
Any scale dependence will introduce structure in the ob- 
servables that, in principle, increases the degeneracies be- 
tween the MGPs and other cosmological parameters such 
as the spectral tilt. However, including high redshift data 
such as CMB measurements in any analysis would break 
the additional degeneracies and we do not expect param- 
eter constraints to be affected significantly if all relevant 
scales are covered by the observations to sufficient accu- 
racy. In this work we will restrict ourselves to elucidat- 
ing the utility, or otherwise, of including cluster counts 
in combined constraints on MGPs. 

The background expansion history is already constrained 
to be close to that of a ACDM model, we will therefore as- 
sume that the modified gravity mimics the expansion his- 
tory of a standard ACDM setup. Having made this choice 
we are then left with probes of inhomogeneity as observ- 
ables that could constrain any modifications. Through- 
out this work we take a fiducial ACDM cosmology de- 
scribed by the following parameters: The dimensionless 
Hubble rate in units of 100 Km s"^ Mpc"\ h = 0.71, the 
density of matter (baryons + dark matter) and dark en- 
ergy in units of the critical energy density, flm = 0.265, 
and r^A = 0.735 respectively, the optical depth to recom- 
bination r = 0.088, the amplitude of primordial, super 
horizon curvature perturbations \og{10^'^ As) ~ 3.071 at 
k = 0.05 h Mpc~^ and their spectral index n^ — 0.963. 
These parameters correspond to the WMAP 7-year best- 
fit parameters [Ml ■ This model yields a large scale struc- 
ture normalisation of ag = 0.804 for the standard devia- 
tions of fiuctuations on scales oi 8h~^ Mpc. 



III. OBSERVABLES 

A number of authors have examined the use of combi- 
nations of observables for forecasting future constraints 
on MGPs [H [T^HT^ . Combinations have included weak 
lensing, CMB, galaxy redshift surveys, peculiar velocity 



surveys and cluster counts. In this work we will exam- 
ine the combination of CMB cross-correlated with weak 
lensing surveys and combined with cluster counts. We 
have made this choice of observables due to the inher- 
ent simplicity in their sensitivity to MGPs and the po- 
tential to unambiguously interpret the data. The CMB 
on largest scales will be sensitive to any modification of 
gravity through the ISW effect as probed by photons, 
whilst on the smaller scales will still provide uncorrelated 
constraints on the conventional parameters of the cos- 
mological model. Weak lensing, being a relatively low- 
redshift sourced signal, will provide direct constraints 
on the MGPs, also, as probed by photons. The cross- 
correlation of the two observables will serve to enhance 
the sensitivity to the MGPs and to reduce the degenera- 
cies between the MGPs and other parameters. 

As with all probes of MGPs, a set of observables that 
depends on the growth of non-relativistic matter pertur- 
bations is required to constrain any difference between 
metric perturbations. For this investigation we have cho- 
sen to focus on cluster counts. These offer the prospect of 
an unbiased tracer of the dark matter distribution as op- 
posed to measurements of galaxy redshift power spectra 
which are known to suffer from scale dependent biasing. 

To calculate the theoretical predictions from models of 
modified gravity we will employ the CAMB sources ^ 
package which has been modified to include modifications 
of gravity during the free streaming regime after recombi- 
nation (MGCAMB [30'). The CAMB sources package allows 
the calculation of full-sky power spectra for a number of 
generic 'sources' which can be those of eg. weak lensing, 
21cm emission, etc. A useful feature of the package is 
the capability of separating out the contribution to an 
individual source into redshift bins, a feature we will be 
making use of for our forecasted weak lensing surveys and 
cluster count predictions. All the models used assume a 
standard ACDM background and the effect of modified 
gravity is solely to change the evolution of perturbations 
around the background. Additionally, We use the mat- 
ter power spectrum calculated by the code as an input 
for our calculation of the mass functions, which are in- 
tegrated to give the expected cluster counts described 
below. 



Throughout this Section and as we discuss our results we 
will refer to a number of 'Stages' of observations which 
summarise a time-line of future surveys in each cate- 
gory of observables discussed. The exact definition of 



A. Cluster Counts 

Galaxy clusters are some of the largest collapsed struc- 
tures in the universe. According to the standard ACDM 
cosmology, they typically consist of hot gas bound in a 
large cold dark matter halo. Clusters have been looked 
at in the context of constraining dark energy [TTJ fSStilU] , 
and some of the studies looking at constraining 7 or the 
MGPs [SlinillHlllIHll]. They are a useful cosmological 
probe as their size corresponds to scales near the linear 
to non-linear transition in the underlying dark matter 
power spectrum. This has several consequences: they 
probe the tail of the matter perturbation spectrum and 
are therefore a sensitive probe of growth. In addition, 
galaxy cluster counts can be predicted accurately from 
linear theory, using semi-analytic formulae or formulae 
calibrated from N-body simulations. These prescriptions 
have been accurately calibrated using ACDM models but 
may need to be revisited if they are to be extended to 
modified dark energy or gravity models. The formulae 
that are calibrated by N-body simulations work over a 
range of cosmologies, but their suitability to perturbed 
dark energy and modified gravity cosmologies have not 
been investigated fully. Some specific dark energy models 
have been looked at gSl US] as well as DGP [13 HI |13- 
I3H] and f{R) [SD]. In most cases the non-linear fitting 
functions and mass functions have been found to be suf- 
ficiently accurate [5T] ^. 

As previously mentioned, we are not taking into account 
any physical screening mechanism in this work which 
would naturally lead to modifications to the mass func- 
tion and any semi-analytical prediction of cluster counts. 
We will look into the impact of screening scales on the 
predictions in future work [^ and use the ACDM cali- 
brated predictions for the forecasting exercise being car- 
ried out here. 

The number of clusters observable over a fraction of the 
sky /sky and with a redshift dependent mass resolution 
limit Miim(z) in a redshift bin spanning the interval z 
to z 4- Az can be calculated by integrating the comoving 
number density dn/dM of objects with mass M 






dn 

Mi„„(^') dM 



dM, (5) 



the Stages considered will be discussed in Section |IV] af- 
ter we introduce our three choices of observables in more 
detail. 



where dV/dzdfl = r'^(z)/ H(z) is the comoving volume at 
redshift z in a flat universe, with H{z) the Hubble rate 
and r{z) — L dz' /H{z') is the comoving distance to that 
redshift. 

Much work has gone into predicting the shape of the mass 
function dn/dM for a given linear power spectrum start- 
ing with the semi-analytical Press-Schechter formalism 
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^ See also [41] for references regarding testing the fitting function 
from [52] 



given by the mass function 
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where p is the background density of dark matter today, 
5c = 1.686 is a critical density contrast and a\i is the 
variance of the dark matter fluctuations in a spheres of 
radius R — (3M/47rp)^/^ defined by the integral of the 
linear matter power spectrum P{k) over wavenumber k 
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with top-hat filter function 



^(^^) - ^ [-jkRT ' ikRr) ■ ^^ 

Successive studies have yielded formulae of increasing 
complexity and accuracy [53-57]. particularly for cluster 
predictions. Other approaches have used N-body cali- 
brated empirical formulae for the mass function. In this 
work we adopt the results of ;52] where the mass function 
is expressed as 
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(9) 

Since we keep an expansion history that is consistent with 
ACDM, the effect of the modified gravity parameters will 
be to change the growth history and hence the matter 
power spectrum. The linear matter power spectrum is 
calculated at the desired redshifts by the MGCAMB code, 
and this is then fed into the mass function ^ and cluster 
abundance ([5]). 

As can be seen in (lo]), the number of clusters is expo- 
nentially sensitive to the amount of growth that has oc- 
curred, and should therefore provide a strong constraint 
on the p parameter. This sensitivity can be seen in 
Figure 111 where we show the expected cluster counts 
in redshift bins of width Az = 0.1 for a Planck-like 
experiment observing 70% of the sky. The curve was 
obtained from J9]) assuming a constant limiting mass 
A^iim = 5.0 X 10 M© , consistent with what is expected 
of Planck [20]. We also show the increase in cluster 
counts expected for the same observational setup with 
p = 1.0023, 1.0046 and 1.0069 corresponding to 1, 2, 
and 3(7 deviations from GR for our Stage I forecasted 
constraints (see Section Vl). These values represent a de- 
parture from the GR value of p of less than 1%, however 
they lead to an increase of between 5 and 10% in the 
number of observed clusters. 

Our theoretical predictions for the number of clusters 
in redshift bins will be compared to predicted SZ cata- 
logues for a number of future observational stages. The 
SZ effect |58] is a nearly redshift independent tracer of 
clusters that is due to the re-scattering of CMB photons 
by hot intracluster gas. The observational limits on SZ 




FIG. 1: Cluster counts in 20 redshift bins for the fiducial 
ACDM model calculated using ([9l with redshift independent 
Ai^iim = 5.0 X IO^^Mq, solid (black). We also show the ex- 
pected increase in counts for models with p = 1.0023 (short- 
dashed, red), 1.0046 (long-dashed, green), and 1.0069 (dash- 
dotted, blue) corresponding to 1, 2, and 3a" deviations from 
GR for our Stage I forecasted constraints (see Section [v|). 



observations are, in principle, determined simply by res- 
olution and sky coverage. The frequency dependence of 
the effect also makes it possible to cleanly identify clus- 
ters in any multi-frequency CMB observations spanning 
the null frequency of ^ 220GHz where the effect changes 
sign. An additional advantage of SZ surveys over X-ray 
surveys of clusters is that the effect is less sensitive to the 
internal structure of the cluster than X-ray emission and 
this should make it easier to obtain unbiased estimates 
of the cluster masses [39l ESj . 

In principle, the limiting mass for an SZ survey is a red- 
shift dependent quantity. The source of this is the change 
of apparent size of the cluster with redshift; the actual 
distortion to the temperature has no redshift dependence. 
We calculated the effect that the full redshift dependent 
limiting mass had on the constraints (see Appendix [A| . 
Since the effect of the redshift dependent limiting mass 
on the constraints turns out to be negligible, we have 
used a constant, redshift independent limiting mass for 
all our forecasts. 

The exponential sensitivity of cluster counts to the ampli- 
tude of the underlying density perturbations introduces 
some issues of accuracy. The mass of a given cluster must 
necessarily be estimated from some proxy signal such as 
X-ray temperature or SZ flux. Use of cluster counts to 
constrain model parameters is therefore subject to any 
bias introduced in the determination of the cluster mass 
from the available information. Any forecasts that do 
not take this uncertainty into account may potentially 




FIG. 2: Total redshift distribution of the sources (solid line) 
and redshift coverage of the four weak lensing bins used for 
forecasting the weak lensing experiments. In calculating the 
redshift coverage of the bins, we have marginalised over the 
uncertainty in determining the redshift of the sources. 



underestimate the errors in model parameters. 
Following |60) . we consider the probability of assigning 

a mass M to a cluster of true mass M to be given by a 
Gaussian distribution in the logarithm of M 



P{M\M)dM = 
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where A^ = InM and A^ = InM. In defining the dis- 
tribution we have assumed there is no systematic bias 
in the estimates of the cluster masses, only a scatter in- 
duced by the uncertainty. The scatter is parametrised by 
the standard deviation of the distribution ct^vi . 
We can now calculate the expected comoving number 
density given the distribution in estimated masses with 
observational cutoff Afii,„ 
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Substituting in the expression for the probability and car- 
rying out the integral yields 
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Integrating the density over redshifts gives a modified 
number count as a function of the uncertainty a^ 



Naz(.cfm) = 47r/sky / dz' 

J z 
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This can be used in our forecasts to account for the un- 
certainty. We will adopt a reference value of aM = 0.25 
[5U] . We will include the a^ as an extra parameter in the 
Fisher matrices for the cluster counts and marginalise it 
out by eliminating its contribution to the inverse Fisher 
matrix. 



B. CMB 

With the release of Planck satellite results only a few 
years away we are entering an era where observations 
of the CMB total intensity spectrum will have reached 
the sample variance limit throughout scales where pri- 
mary effects dominate the signal. As such, in all of our 
forecasts, we will be assuming a Planck-like observation 
of the CMB angular power spectrum out to multipoles 

e = 2000. 

The sensitivity to MGPs in the CMB spectrum is re- 
stricted to the largest scales. This is due the constraint 
that gravity is not modified at early times and through re- 
combination. This means anisotropics generated through 
the Sachs Wolfe effect on super-horizon scales and acous- 
tic effects on sub-horizon scales at recombination will not 
be affected by our late time modifications. The only sig- 
nal will arise on the largest scales due to the ISW effect 
which is sourced as the Universe transitions into a dark 
energy dominated model and the potential starts to de- 
cay. The effect can be described by the integral of the 
time-derivative of the sum of metric potentials along the 
line of sight as photons free stream after recombination 
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where {ST/T)e is the multipole expanded, fourier trans- 
form of the CMB temperature fluctuation at conformal 
time today, tq, is the optical depth, and the spherical 
Bessel functions ji describe the projection of plane-wave 
modes on the celestial sphere. The ISW contributes to 
the power on the largest scales as it is only sourced at 
late times (we will disregard any effect from the early 
ISW effect due to potential evolution close to the radia- 
tion to matter transition). These modes are fundamen- 
tally ill-sampled due to the small number of agm coef- 
ficients on these scales and also suffer degeneracy with 
other effects such as the Sachs Wolfe signal from last 
scattering and the effect of reionisation. One way to 
gain further constraining power from such a sample vari- 
ance limited signal is to cross-correlate the CMB with 
other large scale observables, indeed, low-significance de- 
tections of a dark energy component have been reported 
via cross-correlation with a number of tracer of large 
scale structure [5TH55] . We will also take advantage of 
this by cross-correlating with template weak lensing sur- 
veys. 

As our template CMB observable we take an angular 
power spectrum Ci from the best-fit ACDM model and 



add sample and noise variance according to a Planck-like 
survey covering an area corresponding to 70% of the sky. 
The error at each multipole can be calculated from the 
observational parameters via 



5Ci 



/sky(2^+l) 



{Ct + Ni) 



(15) 



where Ni is a function of resolution (0), the number of de- 
tectors (A'dct) and Noise Equivalent Temperature (Tnet) 
in a given channel, as well as the overall integration time 

(ttot) 



Nf = 
with a evaluated as 



■ exp 



8 In 2 



47r/sky Tnet 



(16) 



(17) 



For our Planck-like survey, we will use the values 6* = 7.1 
arcminutes, iVdet = 32, Tnet = 62.0 ^iK sl/^ /sky = 0.7 
and itot = 14 months. 



C. Weak lensing 

The third observable we will use is the convergence power 
spectrum from weak lensing surveys. Weak lensing is a 
relatively new cosmological tool and is a measure of the 
small distortions of background galaxies due to gravita- 
tional lensing by large scale structure [7D]. Distortions 
of individual background galaxies are virtually impossi- 
ble to measure due to the intrinsic ellipticity of galaxies. 
However, statistical results averaging over large numbers 
of galaxies are now routinely reported |71j . The conver- 
gence K is a measure of the Laplacian of the potentials 
responsible for the lensing along the line of sight and can 
be calculated as an integral over the comoving radial dis- 
tance r as 



$(1 + T])g{r)dr , 



where we have expressed quantities in the Fourier domain 
and g{r) is a filter function determined by the redshift 
distribution of background galaxies being lensed u'[r(z)] 



l{r) = r 



1 



w{r ) dr . 



(19) 



The convergence power spectra can be calculated using 
MGCAMB for the modified gravity case and can be split 
into contributions from separate redshift bins assuming 
the observations are able to obtain sufficiently accurate 
photometric redshifts of the background sources. In all 
cases we will include multipoles i < 2000 to avoid com- 
plications that arise due to lensing from non-linear scales. 



TABLE I: Parameters used for the three Stages of future ob- 
servations used in our forecasts. 





r-K 
Isky 


Ns 


rSZ 
Isky 


Miim (Mo) 


Stage I 


0.121 


2.14 X 10* 


0.7 


5.0 X 10" 


Stage II 


0.485 


3.6 X 10* 


0.7 


2.5 X 10" 


Stage III 


0.485 


2.88 X 10^ 


0.7 


1.0 X 10" 



For our initial weak lensing survey, we consider a DES- 
like survey. DES is a ground based survey at the Cerro- 
Tololo Inter- American Observatory in Chile that is sched- 
uled to begin observations in 2011. It will survey 5000 sq 
deg over 5 years and aims to constrain dark energy with 
4 probes: supernovae, BAO, galaxy clusters and weak 
lensing, the latter being the probe we are interested in 
here. We consider 4 redshift bins between z — and 
z = 2, following the prescription in [12J. 
One of the biggest sources of error in these surveys will 
be errors in the photometric determination of redshifts of 
the background galaxies. We model the overall redshift 
distribution of the sources as 



4 z2 

W{z) = Ng —r= -jr exp 



(20) 



such that Jw{z)dz = Ng, the total number of back- 
ground galaxies and z^, defines the median redshift of the 
distribution. We take a reference value of z* — 0.46 for 
our template weak lensing surveys. 

We can take into account the uncertainty in photometric 
redshifts when breaking down the signal into contribu- 
tions from different redshift bins. We define four, over- 
lapping distributions Wi{z) with i =1, 2, 3, and 4 with 
the constraint 



w{z) =J2wiiz), 



(21) 



i=l 



with 



(18) Wi{z) = -w{z) 



erfc 



Zi-l 



V2a{z) 



erfc 



V2a{z) 



(22) 

The bins are centred at redshifts z^ corresponding to 0.1, 
0.5, 0.9, and 1.3, as shown in Figure [2J with the photo- 
metric redshift error given by u = 0.05(1 -I- z). 
We can model the statistical error in the angular power 
spectrum of the convergence from each redshift bin as a 
sum of sample and noise contributions 



6C7 



/sky(2^+l) 



C^ 



h' 



N'i. 



(23) 



where (7^) is the variance of the intrinsic ellipticity of a 
typical galaxy in the survey and N^ = J Wi{z) dz is the 
number of galaxies in each redshift bin. See table IT] for 



^ 0,0001 ? 




FIG. 3: The cross-correlation signal between CMB temper- 
ature and weak lensing convergence from Stage I type sur- 
vey combination (see Table |l| as a function of angular multi- 
pole I. The curves shown are for the fiducial ACDM model. 
The cross-correlation is shown in four redshift bins with am- 
plitude increasing with redshift. The inset shows the frac- 
tional change in the cross-correlation signal with MGP values 
r\ = 1.03 and ^ = 1.008 corresponding to a Icr shift given our 
final constraints (see Section [v|. 



the parameters used to model the three Stages of weak 
lensing experiments. In all cases (7^) = 0.16 and the 
photometric redshift error is given by u = 0.05(1 -I- z). 
With the advent of large scale weak lensing surveys 
the possibility of cross-correlating CMB and convergence 
maps will become a reality. This will make use of all the 
available information in the data since the signal in two 
such maps will be correlated. The cross-correlation will 
be most useful in this case since dependence on MGPs 
is expected to be strongest on the largest scales. In Fig- 
ure [3] we show the expected cross-correlation signal be- 
tween CMB and the four weak lensing survey redshift 
bins used in our analysis. The inset shows the fractional 
change in the cross-correlation induced by a change in 
MGPs from their fiducial GR values to 77 = 1.03 and 
/i — 1.008. These correspond to la deviations from the 
fiducial values given our forecasted constraint obtained 
from CMB and weak lensing data alone (see Section [V| . 



IV. FORECASTS 

In this Section we carry out forecasts for a number of 
future observational 'Stages'. Since the Planck satellite 
will provide a sample variance limited map of CMB total 
intensity anisotropics covering angular scales where the 
signal of interest lies we will use our Planck-like setup as 
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FIG. 4: Fisher constraints on 77 and /i from the combination of 
Stage I CMB and weak lensing observations with and without 
cross-correlation (orange/dot-dashed ellipse and red/solid el- 
lipse respectively). 



the CMB contribution throughout. For weak lensing and 
cluster counts we will assume three distinct observational 
stages corresponding to short, medium, and long-term 
development of survey sizes and accuracies. 

• Stage I: Corresponds to a Planck-like cluster sur- 
vey and a DES (Dark Energy Survey)-like weak 
lensing survey. DES is being carried out on the 
Cerro Tololo Inter-American Observatory in the 
Chilean Andes and should start taking data in late 
2011. DES is modelled with a total of 214 million 
galaxies over 12% of the sky. 

• Stage II: This stage includes a weak lensing sur- 
vey based on PANSTARRS and as such is modelled 
with a total of 360 million galaxies over 48.5% of 
the sky. The SZ survey is modelled by keeping the 
same /sky as with Planck, but lowering the limiting 
mass to 2.5 x lO^^M©. This corresponds to lower- 
ing the smallest change in flux that the SZ survey 
can detect. 

• Stage III: This includes a weak lensing survey 
based on the LSST, due to begin taking data in 
2020. This survey is modelled with the same frac- 
tion of the sky as PANSTARRS (48.5%), but with 
a total number of galaxies of 2880 million. The 
third stage SZ survey assumes a limiting mass of 
1.0 X IO^Mq. 

The sampling characteristics for the three Stages are 
summarised in Table IT] We have kept the same intrin- 
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0,96 0,98 1,00 1.02 1,04 



V 



FIG. 5: Fisher constraints on rj and jj, from cross-correlated 
CMB and weak lensing measurements (Stage I) are shown 
by the red area (dot-dashed contour). The improvement ob- 
tained by adding Stage I cluster counts is seen in the green 
area (solid contour). The outer ellipse corresponds to the 
inner ellipse of Fig. [4] Adding an uncertainty in the mass as- 
signment for the clusters of aM = 0.25 decreases the impact 
of adding cluster data as shown by the yellow (dashed) ellipse. 



sic ellipticity, redshift bins and photometric errors as- 
sumptions for all Stages as it is unclear how photometric 
redshift resolution will evolve as surveys increase in size 
and complexity. 

Our forecasts are based on Fisher matrix estimates of er- 
rors in a subset of parameters that comprises the MGPs, 
77 and /i, and the two parameters from the standard model 
that are expected to be most correlated with them, the 
total matter density ilm and the primordial amplitude of 
scalar curvature perturbations A. In varying only these 
four parameters we are assuming that the remaining four 
standard ACDM parameters are well constrained by sig- 
nals orthogonal to those being used in our analysis. A 
combination of high resolution CMB spectra, including 
polarisation i?- modes, and a standard prior on the value 
of the Hubble rate will ensure most of the remaining pa- 
rameters are fixed to within a few percent of their nom- 
inal values, which should have a minimal impact on the 
four parameters being considered in this work. 
An alternative to our Fisher matrix based method would 
be to Monte Carlo Markov Chain (MCMC) sample the 
joint posterior distribution of our parameter set by eval- 
uating the likelihood explicitly. However this would in- 
volve the definition of a likelihood as a function of real- 
isations of the observables (and their cross-correlation) 
and the added complexity is not warranted for this kind 
of exercise at this stage. 



It should be noted that use of a Fisher matrix to estimate 
parameter errors assumes that the observables are dis- 
tributed as Gaussian variates. This is not true in all cases 
considered here since we are considering power spectra 
and number counts but as long as the true answer lies 
close to our fiducial values for the four parameters the 
errors should give a good indication of the constraints. 
Given a set of n, uncorrelated measurements C'a, with 
a = 1, ..., n and measurement errors SCa, the Fisher ma- 
trix for a set of m parameters Aq with a = l,...,m, can 
be evaluated as 



FaP — 



^ 1 dCgdCg 1 

^ 5Ca dXa dXp SCa 



(24) 



where /? = 1, ...,m. 

The Fisher matrix represents the ensemble average of the 
negative curvature in the log likelihood of the model pa- 
rameters and its inverse, in this limit, is therefore the 
covariance matrix in those parameters. The Fisher ma- 
trix is simple to evaluate since it involves only the first 
derivatives of the signal with respect to the model pa- 
rameters. These can be evaluated either analytically or 
numerically. Here we use a central difference scheme to 
numerically approximate the derivatives to second order 
in the step-size. The central difference is sampled by 
evaluating the models with given step-sizes either side of 
the fiducial model in all parameter directions. 
For the cluster counts case the measurement consists of 
counts in each of twenty redshift bins Ni, as shown in 
Figure [l] for the fiducial model, and the Fisher matrix is 
calculated as 



^ap 



1 dN, dN, 



-^ cr^^ dXa dXp 



E 



(25) 



where we have assumed a shot noise model for the error 
in the counts. 

Fisher matrices from independent data can be added and 
then inverted to obtain error estimates for the combina- 
tion of data. Since we will rely heavily on the cross- 
correlation of CMB and weak lensing measurements to 
extract the relevant signal we calculate the combined 
Fisher matrix for this cross-correlation to add to ( 25 1 . 



The Fisher matrix formalism can be easily extended to 
the correlated measurement case. We treat the combina- 
tion of CMB and weak lensing measurements at a given 
multipole ^ as a matrix, Cg, of angular, cross-correlation 
power spectra with dimension n x n. In our case the in- 
dex n spans both CMB and convergence angular power 
spectrum measurements (over the four redshift bins) i.e. 
T, Ki, K2, K3, Ki with the symmetric form 



fC^> 



■TT 



c: 



Tki 



c. 



'Tk2 



\ 






(26) 
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Stage I 


Stage II 


Stage III 


Parameter 


CMB ® WL 


CMB ® WL + CC 


CMB ® WL 


CMB ® WL + CC 


CMB (g) WL 


CMB ® WL + CC 


A-q 
Am 


3.0 X 10-2 
8.4 X 10'^ 


5.5 X 10-^(1.5 X 10-2) 
2.3 X 10^^(4.5 X 10"^) 


1.8 X 10-2 
5.5 X 10-=* 


4.4 X 10-^(8.8 X 10-^) 
2.1 X 10^^(3.1 X 10"^) 


1.5 X 10-2 
4.8 X 10"^ 


2.4 X 10-^(5.5 X 10-^) 

1.5 X 10^^(2.1 X 10"^) 



TABLE II: Fisher errors on rj and fi for Stage I, II, and III observations for different combinations of observables: cross- 
correlation of CMB and weak lensing (CMB (g) WL) and addition cluster counts (CC). The values in brackets show liow the 
errors are affected by the marginalisation over a mass assignment error. 



Each measurement matrix will have a corresponding co- 
variance matrix with non-zero off-diagonal contributions 
due to the correlations in the signal part of the measure- 
ments whilst the diagonals will have contributions from 
both signal and noise 



f{5C, 



S,= 



((5Cf''^)2 



V 



\ 



/ 



(27) 



where the diagonal elements correspond to the square of 
standard deviations (15 1 and (23) and the off-diagonal 



terms can be evaluated using 



{5C, 



XY\2 



(2^+1) 



r<xxr<yy 



fX fY 

•I sky -' sky 



{c, 



XY\2 



min(/sky./sky 



(28) 

Here Cf^ is the model cross-correlation power spectrum 
for the two observables which is also computed by CAMB 
sources. 

The Fisher matrix for this generalised case can be eval- 
uated as 



F, 



ap 



ETr 



da 



•s; 



dp 



(29) 



A further contribution to the Fisher matrix is given by 
the sample variance depending on the parameters 



F, 



all 



jETr 



da 



. 1 as, ^_i 



9/3 



(30) 



This term is sub-dominant to the first, however we have 
included it in the analysis for the sake of completeness. 

V. RESULTS 

We initially considered the case where the MGPs switch 
from their fiducial, GR value to a new, time independent 



value, at a fixed redshift. 



There is a wide range 



in redshifts, spanning from the time of recombination to 
the end of the so called 'dark-ages', where modifications 
to GR could come into effect. Here we chose z^g = 30 
for the time-independent case which ensures that most 
non-linear structures observed formed in the presence of 



f^ l.OQO--- 




0.994^ 
0.98 



FIG. 6; Fisher constraints on 77 and fj, from the combination 
of CMB, weak lensing, and cluster counts for the three ob- 
servational stages. Stage I, red (dot-dashed); Stage II, yellow 
(dashed); Stage III, green (solid). 



the modifications. We also consider the case where the 
MGPs have a simple, linear redshift dependence between 
their GR values at Zmg and their late time values rjo and 
fiQ at z — 0. Specifically we show how our results change 
for a choice Zmg = 3 and Zmg = 1- 

The CMB is most effective in constraining the standard 
parameters of the concordance model of cosmology. The 
only effect of the MGPs is on the ISW effect, and this 
has several shortcomings. The increase of power due to 
the MGPs is not large, and this is compounded by the 
higher cosmic variance on these scales. The uncertain- 
ties in weak lensing mostly come down to a degeneracy 
between i] and A. This, is why, when combining CMB 
and Weak lensing observations, the constraints on the 
modified gravity parameters improve significantly. The 
extra information obtained from cross correlating the two 
observations emphasises this complementarity of the two 
sets of data. Figure |4] shows how the constraints on rj 
and /I improve when the cross-correlation between CMB 
and weak lensing is included in the Fisher estimates. The 
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FIG. 7: Same as in Fig. IS] but with a mass assignment uncer- 
tainty of (jm = 0.25 included in the constraints. 



contours shown the 1-cr constraints obtained from the in- 
verse of the Fisher matrix for all four parameters. We 
ensure our final errors account for marginalisation over 
our nuisance parameters Vim and In A by only ignoring 
the matrix elements corresponding to these parameters 
after inverting the full Fisher matrix. 
Even after taking into account the cross-correlation be- 
tween CMB and weak lensing a significant degeneracy in 
•q remains due to the correlation between the two gravita- 
tional potentials in both signals. Adding cluster counts 
to the mix improves the situation since their signal is 
highly sensitive to just ^ and therefore to the parame- 
ter /i, which affects the growth rate of structure. This 
also has a large effect on the r] parameter as the addition 
breaks the remaining correlations and allows the CMB 
and weak lensing observations to constrain 77. The ef- 
fect of adding Stage I cluster counts is shown in Figure [5] 
As mentioned earlier, there are some uncertainties asso- 
ciated with cluster counts. As can be seen readily from 
Figure |5l although marginalising over these uncertain- 
ties reduces the impact of clusters, they still add to the 
constraining power of the CMB and Weak lensing. 
The comparison for Stage I, II, and III Fisher results in 
77 and /x are shown in Figure |6] The effect of adding a 
marginalisation over the mass assignment uncertainty is 
shown in Figure [7] The constraints worsen by a factor 
of between 2 and 3 when the uncertainty is taken into 
account but we stress that our estimates are conservative 
since we have not allowed for any improvement in um in 
successive Stages. 

A summary of the results is shown in TablelTTlfor all three 
Stages. The forecasted errors A?7 and A/^ are shown with 



FIG. 8: Fisher constraints on 'q and ^ from the com- 
bined CMB, weak lensing, and cluster counts for the 
time-independent (solid/yellow) and time-dependent (dot- 
dashed/red) modified gravity parameters. In both cases 
^mg = 30 but for the latter a linear redshift-dependence is 



assumed (311 



and without addition of cluster counts (CC). All con- 
straints use the cross-correlation signal between CMB 
and weak lensing (CMB ® WL) and the values in brack- 
ets show the effect of adding in the mass assignment un- 
certainty to the cluster counts. 

In general, we find that the MGPs can be constrained to 
within a few percent using our forecasted surveys. The /i 
parameters is most constrained by the observations. This 
is not surprising since it is the parameter that affects the 
growth history directly. Including the mass assignment 
error in the cluster counts reduces the constraints by ap- 
proximately a factor of 2. The improvement obtained 
from successive Stages is limited, however the Stages re- 
cover the constraining power that was lost to the mass 
assignment error of the cluster counts. This reinforces 
the need for increased sensitivity and survey size for fu- 
ture observations as it hedges against systematic effects 
such as mass and photometric redshift resolution taken 
into account here. 

To gauge the robustness of our predictions with respect 
to our choice of fixed 77 and /x, we add a linear time de- 
pendence in the MGPs of the form 



/i(z) = (/lo - 1) 



(31) 



with a similar expression for 77. Here, /xq or 770 is the 
value of the MGP today (z = 0), which is the parameter 
that will be constrained in the Fisher matrix analysis. 
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As such, the choice of Zmg, which is not included as a 
parameter in the Fisher analysis, fixes the gradient of 
the linear time dependence. We are therefore testing our 
results with respect to a choice of Zmg rather than ex- 
tending the parameter space to include a parametrised 
time dependence. 

As before, z^g is the redshift at which the modification 
to gravity 'switches on'. As shown in Figure [8] allow- 
ing for this simple time dependence in the MGPs does 
not have a large impact on constraints for the case where 
■Zmg = 30. Although the time dependence reduces the 
impact of the MGPs on the observations for the same 
value of the MGPs today, the relative effect is also differ- 
ent in each of the signals and leads to a small change in 
degeneracies. An additional signal is sourced in this case 
since the time dependence of rj will affect the late-time 
ISW. These combination of effects partly cancel, which 
results in a small effect on the final constraints. It should 
be stressed that there are many different choices of time 
dependence for the MGPs and (31) is only one of these. 



In principle some other form of time dependence could 
lead to much larger effects in this type of comparison. 
A theoretically motivated for for the time dependence 
would be better justified, however, this would build in 
a model dependence which we are explicitly avoiding in 
this work. 

Our results for Zmg — 30, 3, and 1 are summarised in 
Table |III| w hich can be compared to the Stage I results 
of Table |II| We find the effect becomes large only at rel- 
atively low redshifts with the choice of z^g = 1 leading 
to an order of magnitude increase in Fisher errors. This 
result is not indicative of the ability of the observations 
to constrain any possible time behaviour of the modified 
gravity theory but rather it illustrates the sensitivity of 
the observables to the departure from GR integrated over 
time. 



VI. DISCUSSION 

Future generations of large area surveys hold much 
promise in testing the validity of GR. In this work we 
have examined the impact, on MGP constraints, of fu- 
ture cosmic shear and cluster count data used in com- 
bination with CMB measurement. We have found that 
this combination can constrain the MGPs to sub-percent 
accuracy. In particular the inclusion of cluster counts, 
which are highly sensitive to any change in the growth 
of the matter perturbations, adds a strong refinement in 
the search for any deviations from the standard GR val- 
ues of the MGPs. However we have also shown that the 
inclusion of a simple model for calibration uncertainties 
in the counts can affect the errors. 

Our Fisher matrix analysis has shown that the MGPs 
may be statistically constrained to fractional levels of be- 
low a percent. This clearly supercedes the current mod- 
eling accuracy of weak lensing and cluster abundances 
in particular. In the case of weak lensing, even at large 



scales, non-linear contributions to the matter power spec- 
trum are not modeled accurately enough whilst in the 
cluster case the mass function formalism is certainly not 
accurate at the 1% level. Any analysis of actual data 
will necessarily have to include a marginalisation over 
these and other modeling uncertainties. This can be done 
either by including the model uncertainties in the final 
Fisher matrix and marginalising by cutting out the di- 
rections in its inverse or by including the model parame- 
ters as extra directions in MCMC methods. Future MGP 
constraints may well be limited by such systematic effects 
rather than the statistical limits of the observations. 
There are, of course, other observables that can constrain 
the growth of matter perturbations such as large scale 
structure, peculiar velocities, redshift space distortions 
and future high redshift 21cm surveys. These could all 
replace our choice of cluster counts as a signal used to 
break the degeneracy between gravitational potentials. 
All of these, however, present a number of problems ei- 
ther in their interpretation as biased tracers or in the 
technological challenges involved in large scale surveys. 
The observations of clusters on the other hand is a rel- 
atively simple procedure and SZ surveys of clusters are 
a necessary by-product of current and future CMB ex- 
periments. The only difficulty involved in this signal is 
the theoretical modelling required to use them to obtain 
constraints due to the uncertainty in the choice of mass 
functions to be used for theories of modified gravity. 
In obtaining our forecasted constraints we have only em- 
ployed the simplest description of phenomenologically 
modified gravity. The model building for such theories 
has already shown that there there is much more com- 
plexity to explore. All models must incorporate a screen- 
ing mechanism to satisfy standard solar system and lab- 
oratory tests of gravity. Screening mechanism are only 
one way in which modified gravity models may include 
significant scale dependence, which we have not taken 
into account in this work and may provide more testable 
predictions and therefore more ability to differentiate be- 
tween models. Concrete models may also lead to mod- 
ifications to the background expansion, which we have 
not explored here. This can also lead to additional con- 
straints on model parameters. 

Allowing for a very general scale dependence of the MGPs 
will inevitably introduce a scale dependence into the con- 
straints. In particular our combination of data generates 
constraints based on large scales, due to the ISW effect, 
and a fairly narrow range of scales of around 10 Mpc 
from cluster counts. If there is a significant scale de- 
pendence of the MGPs between the two scales then the 
combination of data considered here would give stronger 
constraints on some scales than others. A complete treat- 
ment would include a Fisher eigenmode analysis for an 
extended range of parameters describing the scale depen- 
dence of the MGPs [H m [H]. We leave this for future 
work for the combination of data considered here. 
The experimental outlook is promising for the observa- 
tions we have dealt with in this work. Over the next 5- 



13 



Stage I 
Parameter 



CMB®WL| CMB (g) WL + CC CMB (g) WL 



2.6 X 10"^ 6.4 X 10"^(1.5 x 10"^) 3.3 x 10"^ 
9.1x10"^ 3.5x10^^(5.6x10"^) 1.0x10"^ 



= 3 

CMB (g) WL + CC 



1.9 X 10"^(2.6 X 10"^; 
5.0 X 10^^(7.6 X 10"^; 



CMB (g) WL 



CMB ® WL + CC 



1.1 X 10^^(1.1 X 10"^) 
5.1 X 10^^(5.1 X 10"^) 



A-q 
Afi 



1.1 X 10" 
5.1 X 10" 



TABLE III: Fisher constraints for time-dependent MGPs for Stage I observations. 



10 years, deviations from GR should be well constrained, 
and the concordance cosmology will either be more se- 
cure or may even have undergone a paradigm shift. A 
null result would support the concordance cosmology, a 
conclusion that would be even stronger if dark matter 
had been detected non-gravitationally by then. A detec- 
tion of a deviation from GR would be potentially more 
interesting but would require a completely new theoret- 
ical framework and trigger a search for an underlying 
model for the modifications. Of course, a dark energy 
model with perturbations may turn out to fit the data 
just as well and it is not clear at this point whether this 
degeneracy will ever be broken by observations. 
Acknowledgements We thank Filipe Abdalla for pointing 
out the work on self calibration of clusters to us. This 
work was supported by an STFC studentship. 



Appendix A: Redshift dependent limiting mass for 
cluster surveys 

The full treatment of the limiting mass here is based on 
[75] . The limiting mass at redshift z, for a detector at 



dimensionless frequency x with a flux limit of 5'iini is given 
by 



Mr, 



lim 



Slinid\{z) 



fix)fbAszil + z) 



f^-L 



(Al) 



Here, (iyi(z) is the angular diameter distance, /h is the 
baryon fraction and /(x) is the frequency dependence of 
the SZ effect, given by 



/(^) 



(e^ - 1)2 



a;(e^ + l) 



-4 



ksTc 



The dimensionless frequency is given by x = ^ ^ 
with TcMB the temperature of the CMB. For our analysis 
based on Planck, we used S'ljm = 30 mJy and a frequency 
V — 353 GHz. Both Asz and /3sz are parameters from 
the mass-temperature relation for clusters. The values 
of the constants used here are as follows: /icm — 0.06, 
/3sz = 1.75, and Agz = 3.781 x 10^. 
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